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Abstract 

We formulate a method to solve the coordinate space Faddeev equa- 
tions for positive energies. The method employs hyperspherical coordi- 
nates and analytical expressions for the effective potentials at large dis- 
tances. Realistic computations of the parameters of the resonances and 
the strength functions are carried out for the Borromean halo nucleus ^He 
(n+n+Q) for J" = 0±,1±,2±. 

PACS numbers: 21.45. -fv, 11.80. Jy 31.15.Ja, 21.60. Gx 

Introduction. The three-body continuum problem has been the subject of 
numerous investigations Tremendous progress has been achieved, but still a 
number of problems remain . Many approximate solutions have been invented 
without an emerging established general procedure. Different treatments are 
usually needed for short-range and long-range interactions and for energies below 
or above possible two-body thresholds |^, D ■ It is necessary, but not always 
easy, to distinguish between inaccurate numerical results and shortcomings of 
the basic interactions. 

During the last decade a new class of weakly bound three-body systems, 
nuclear halos, attracted enormous attention 0, §]• If no binary subsystem 
is bound, they are called Borromean nuclei. These concepts are general and 
of interest in many subfields of physics ^, |l^. Accumulating data from such 
systems demand analyses heavily relying on the properties of their continuum 
spectra [0 |l^. However, technical difficulties related to the precise behavior 
at large distance are substantial and so far unsolved. 

Recently a new method with explicit analytical treatment of the large dis- 
tances jl^, was used to obtain bound-state solutions to the Faddeev equa- 
tions. The method is very powerful as seen by the successful investigation of 
the Efimov effect g ||. The purpose of this letter is to generalize the method 
to obtain continuum state solutions. In order to illustrate the efficiency of the 
method we perform a realistic computation of a three-body Borromean halo 
system. 
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Method. The fc'th particle has mass irik and coordmate r^. The two-body 
potentials are Vij. We shall use the three sets of Jacobi coordinates (xi,yi) 
and the corresponding three sets of hyperspherical coordinates (p, at, ^xi, ^yi) 
H 0- The kinetic energy operator is then 

T^—(~n^^/^ — n^^^ + — + —\ (1) 




A2 = — ^|^sin(2a) + — ^ + — ^-4, (2) 
sin(2a) oa'' sm a cos"^ a 

where the angular momentum operators l'^ and ly are related to the x and y 
degrees of freedom. 

The total wave function is now expanded in a complete set of hyperangular 
functions 

po/^ sm(2Q!j 

n=l i=l ^ ^' 

Each of the three components (j)n^ is expressed in the corresponding system of 
Jacobi coordinates and they satisfy for each p the three Faddeev equations 

(A2 _ A„) + ^p'V,, y = 0, (4) 

^ "^sin(2a0 h^^ ^^^sin(2a/) ^' 

where j, fc} is a permutation of {1, 2, 3}. In the absence of bound subsystems 
the eigenvalues A„ approach at large distances the hyperspherical spectrum 
obtained for Vjk = 0, i.e. A„(p oo) = Kn{Kn + 4), where Kn is odd or even 
natural numbers depending on the parity. 

The expansion coefficients fn{p) satisfy the equations 

/ 92 A„ + 15/4 2mE\ 

[rW^^^ (5) 



n' 



The coupling terms P and Q approach zero at least as fast as p ^. For 
Borromean systems we can then choose those solutions ^n' to eq. where the 
large-distance [p 00) boundary conditions for ^ are given by 

/("') (p) -> S„^„,Fi-\Kp) - S,,,,>fW (^p) , (6) 
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where 
order by 



2 f i) 

Zm E/n-" and F^i' are related to the Hankel functions of integer 




2'KKh' 



m 



■2exp ±iKp± —{Kn + -) 



(7) 



The continuum wave functions are orthogonal and normalized to delta functions 
in energy. 

By diagonalization of the S'-matrix we obtain eigenfunctions and eigen- 
phases. The phase shifts reveal the continuum structure of the system. In 
particular, a rapid variation with energy indicates a resonance. A precise com- 
putation of resonances and related widths can be done by use of the complex 
energy method, where eq.(H) is solved ior E = Er ~ iV /2 with the boundary 

condition — 5n.n' \J'^^ H^k\2{'^p) ■ These solutions correspond to poles of 
the S'-matrix Q. 

Large-distance behavior. Eq.(H) can be solved for large distances, where 
for short-range potentials all partial waves, except s-waves, decouple. We ex- 
pand each component on the hyperspherical basis with the quantum numbers 
{lx^ ly, L, Sx, Sy, S, J} where L,S and J are the total orbital angular momentum, 
total spin and total angular momentum, respectively. We express two of the 
Faddeev components (j, k) in the coordinates related to the third Jacobi set (?) 
and project out the partial wave with a given set of angular momentum quantum 
numbers. This operation, leading from the i'th to the j'th Jacobi coordinates, 
is denoted by Rij. 

For large p only small a contribute to the terms proportional to Vjkiri) in 
eq.(^). This is due to the assumption of short-range potentials and because 
Ti oc psinaj. Let us first exphcitly consider the three coupled components, 
characterized by Ixi = and lyi = L and therefore with the same total 
orbital angular momentum L and furthermore with the same spin structure. 
We expand in powers of ai and find the leading order contribution from the 
transformation of such terms to be 



where p is the parity of the permutation k}. Non-zero /^i-values had pro- 
duced higher powers of ai in eq.(0). Thus, the eigenvalues A related to the other 




sin(2(pj_i) 



(8) 
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partial waves decouple at large distances and approach the hyperspherical spec- 
trum. These waves assume the asymptotic behavior on a distance scale defined 
by the range of the interactions. On the other hand the s-waves couple and feel 
consequently the interactions over a distance defined by the scattering lengths. 

We shall now concentrate on a system consisting of two neutrons and a spin- 
zero core. This model directly applies to ^He, a halo nucleus for which a large 
amount of experimental data exists. The model is also a good approximation 
for another halo nucleus, ^^Li [||. 

Due to the antisymmetry between neutrons the three coupled components 
{Ixi = 0, lyi — L, i = 1,2,3) reduce to two and the angular Faddeev equations 
eq.(Q) are to leading order in a (large p) given by 



L{L + 1) 
da\ cos^ ai 



+ p v^^ipsinai) - v (10) 



X(j}-^\p,ai) = -2ai{-l)^ p^vn^{p sin ai)&p , 

+ P^WNc(psinQ!2) - v"^ ) (11) 



9^ ^ L{L + l 
da\ cos^ a2 



2 

L „2„ 



X(/)}^^(/9, Q;2) = -2a2(-l) p WNc(psinQ!2)C^ , 

where ^ A + 4, t'NN(a;i) = V23(a;i/p23)2TO/?i^, V'^c{x2) = Vi3{x2 / p-i'i)'^m / Ti' , 
^ sin(2^) ' ^ sin(2(p) sin(2^) ' ^ ' 

with ip = (pi2, ip = (y523- 

For large p the short range potentials p'^v{psmai) vanish for all ai except in 
a narrow region around zero. Due to this rescaling the effective range approx- 
imation becomes better with p increasing and therefore any potential with the 
same scattering length and effective range would lead to the same results. Let us 
then in the region of large p use square well potentials Vjk{r) = —V^^''Q{r < Ri), 

or equivalently expressed by the reduced quantities Vjk{x) = —Vq''Q{x < Xi = 
RiPjk), where the parameters are adjusted to reproduce the given two-body 
scattering lengths and effective ranges. The corresponding solutions are then 
accurate approximations to our original problem at distances larger than 2Ri 

[0. 

The potentials v{psm.ai) are zero when ai > a^"^ = a.TCsm{Xi/ p). Then 



eqs.(|lO[]ll|) are especially simple, i.e. 
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and the solutions, vanishing at ai — n/2, are given by 



Pl{v, a) = cos a 



lJQ 1 



\ da cos a ^ 











1 sm 









(14) 
(15) 



The potentials w(psinai) are for large p only finite when a,; < ctg*^ ^ 1. 
Then eqs.(|l^) and (0) are 



da 



(16) 



where the wave functions in in eq. (|l|) must be 
eq.(p^) are then 

4^^^)(p,a) = Bf sin(«:,a) - 2a(-l)^^^ci*) 
«2 = _[L(L + 1 



The solutions to 



2 (i) 

p vy -v 



L ' 

,2l 



(17) 
(18) 



Matching the solutions, eqs.(|lj) and (|l^), and their derivatives at ai = a, 
gives a linear set of equations for A^-* and B^^ . Physical solutions are then only 
obtained when the corresponding determinant is zero. This is the quantization 
condition for A and the eigenvalue equation determining the asymptotic behavior 
of A(p). 



Realistic computations for ^He. The practical implementation of the me- 
thod is tested on ^He considered as two neutrons and a ^He-core. The two- 
body interactions reproduce accurately the s-, p- and d-phase shifts up to 20 
MeV. Furthermore, a diagonal three-body force, exp(— /9^/&§), is added in 
eq.(|^) for fine tuning. The range of the three-body force is by its definition 
given in terms of the hyperradius. For ^He, p—2 fm and 3 fm correspond 
roughly to configurations where the neutrons respectively are at the surface 
of the a— particle and outside the surface by an amount equal to their own 
radius. The idea of using the three-body force is to include effects beyond those 
accounted for by the two-body interactions. 

Several phase equivalent parametrizations are possible for each radial shape 
of the two-body potential. They differ in the number of two-body bound states 
of which the lowest s-state is occupied by the core neutrons and therefore subse- 
quently has to be excluded in the computation. The results are very close after 
fine tuning by use of the three-body interaction . We shall therefore only use 
the potentials without bound states. All possible s-, p- and d- waves are included 
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Figure 1: The lowest angular eigenvalues A„ as function of p for angular momentum 
= 1~ (solid lines), and 2^ (dashed lines) for ^He. The dotted lines are the large- 
distance asymptotic behavior. The neutron-neutron and the neutron-^He interactions 
are from @ with (Vi'\rJi'') = (48.2 MeV, 2.33 fm), (-47.40 MeV, 2.30 fm), (-21.93 
MeV, 2.03 fm) for s,p and d waves respectively. The spin orbit parameters are {Vso,rso) 
= (-25.49 MeV, 1.72 fm). Maximum A"„-values up to 142 are used. 



whereas other waves can be ignored to the accuracy we need. The number of 
Jacobi polynomials in the basis expansion is carefully chosen to give accurate 
numerical results up to a distance, typically around 40 fm, where the asymp- 
totic behavior is reached and from then on the asymptotic solutions eqs.(p^ 
and eq.([l^) are used. 

The accurate low-energy continuum spectrum calculations require integra- 
tion of the radial equations up to distances of the order of ten times the sum 
of the scattering lengths. For the n+n-|-a system this is about 180 fm. Too 
small basis size and too small maximum distance are both disastrous for the 
numerical reliability. 

In Fig. 1 we show the two (strictly decoupled ) angular eigenvalue spectra 
for 1~ and 2+. The structure is complicated at small distances where avoided 
level crossings are seen. The lowest level has in both cases an attractive pocket 
unable to bind the system, but still responsible for several resonances. At large 
distance the structure is simpler as the hyperspherical spectrum is approached. 
In the computation we use the asymptotic behavior, also shown on Fig.l. This 
improvement of the procedure is absolutely essential when accurate results are 
required. 

The phase shifts for the cases in Fig. 1 are shown in Fig. 2. The rapid 
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Figure 2: The eigenphases corresponding to the lowest A-values obtained after di- 

agonalization of the 5- matrix for — \~ (solid lines) and J'^ — 2^ (dashed lines). 
The interactions are the same as in Fig.l where a diagonal three-body interaction, 
ySs exp(— p^/6§) with S3 = —31 MeV, 63 = 2.061 fm, is added in all partial waves. 



variation and the subsequent crossing of 7r/2, seen at four energies, are the 
traditional signs of resonances. On the other hand it is also possible to have 
resonance- like structures without phase shifts crossing 7r/2. They may still show 
up as poles of the S-matrix. 

In Table 1 we give for a few spins and parities the two lowest S'-matrix poles 
obtained by the complex energy method. The radial equations were integrated 
numerically up to PmarE=180 fm where the 1/ tail of the effective potential 
becomes negligible. The numerical solutions were then matched at Pmax with the 
Hankel functions i? . Precisely at the pole only the function must match 
the numerical solution. We show the results for two different three-body forces, 
i.e. fine tuned to the ground state energy and to the 2+-resonance. Although the 
differences appear to be relatively small, they are important for the observable 
properties. The two cases in Table 1 can be considered to give the realistic range 
of the possible variation of the three-body force. For angular momentum 0~, 1^ 
and 2~ the relatively small pocket in the effective radial potential combined 
with strong centrifugal barrier hinders the effect of the three-body force on the 
phase shifts and the resonance properties. Apart from the lowest 2+ state all 
these resonances reside above the effective centrifugal barrier and must therefore 
correspond to rather smooth structures in the cross sections. 

An observable less sensitive to the large-distance behavior is related to the 
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Table 1: The real and imaginary values {Er,r) (in McV) of the two lowest S'- matrix 
poles E = Er ~ ir/2 for ^Hc for various spins and parities. The interactions used are 
the same as in Fig.l. The three-body interaction parameters are 6*3 = —7.55 MeV, 
fes = 2.9 fm and S3 = — 31 MeV, 63 = 2.061 fm respectively for the first two and the 
last two columns. Correspondingly the excitation energies are E* = Er + 0.95 MeV 
and E* =Er + 1.54 MeV. 
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0.37 


1.23 


0.45 


0.845 


0.093 


1.05 


0.40 


2- 


0.90 


0.34 


1.82 


0.57 


0.90 


0.34 


1.82 


0.57 



excitations from the ground state. We show in Fig. 3 the lowest three strength 
functions, Sq+^j„{E), as functions of energy both for plane waves and for the 
proper continuum wave functions. We find 91%, 60% and 70% of the strength 
below 5 MeV, respectively for monopole, dipole and quadrupole excitations. 
We notice the usual rise from zero to a maximum and the fall off towards zero 
at large energy. The peak is very pronounced for 2+ reflecting the observed 
resonance of width 0.11 MeV at 0.82 MeV. Above the smooth plane- wave back- 
ground for 1~ is seen a peak at about 0.95 MeV and a shoulder at about 1.8 
MeV. This significant 1~ enhancement is the result of a combination of two 
overlapping broad resonances, see the S-matrix poles in tabic 1. It should be 
detectable although in the same energy region as the 2+-resonance. The nu- 
clear 0+ strength function resembles the plane wave result reflecting broader 
underlying structures. 

Conclusions. Wc have formulated a method to compute low-energy three- 
body continuum spectra for arbitrary short-range potentials. It is based on a 
recent successful method used to calculate bound states by solving the Fad- 
dccv equations in coordinate space. The angular part of the equations are 
treated purely numerically at short distances, whereas the large-distance be- 
havior of eigenvalues and eigenfunctions is computed essentially analytically for 
all partial waves. Combining the results from these two regions allow accurate 
computations up to very large distances. Realistic computations for ground 
state properties, transition matrix elements, phase shifts, resonance energies 
and widths of J"" = 0±.1±,2± are carried out for the Borromcan halo nucleus 
^He. The established J'^ = 2+ resonance is found together with a number of 
other broader resonances. 
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Figure 3: The strength functions, S'-^^E) = | (fi^''l|A^(^^A)||0+> |^ for «He 
as function of energy for transitions from the ground state to 0^ (dotted), 1~ 
(solid) and 2"'" (dashed) excited continuum states. The operator is M{E\,ii) — 
and gir^VA^(»^i) respectively for A = and 1,2. The units are the cor- 

responding sum rules (0+|p*|0+> - (0+|p2|0+>2=749 fni'' for A = and g^(2A + 
l)(0+|ri^|0+)/(47r)=1.31e^fm^ and 8.19e^fm* respectively for A = 1,2, where = 2e 
is the ■'He-charge and ra is the ''He-distance from the ''He center of mass. The inter- 
actions are the same £is in Fig.l with a diagonal three-body interaction added in all 
partial waves. The parameters are S3 = —7.55 MeV, 63 = 2.9 fm for 0"*" and S3 = —31 
MeV, 63 = 2.601 fm for 1" and 2+. The smooth curves correspond to plane waves for 
the continuum states. 
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